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NODAL LENGTH OF STEKLOV EIGENFUNCTIONS ON 
REAL-ANALYTIC RIEMANNIAN SURFACES 

IOSIF POLTEROVICH, DAVID A. SHER, AND JOHN A. TOTH 


Abstract. We prove sharp upper and lower bounds for the nodal length of Steklov eigen¬ 
functions on real-analytic Riemannian surfaces with boundary. The argument involves fre¬ 
quency function methods for harmonic functions in the interior of the surface as well as the 
construction of exponentially accurate approximations for the Steklov eigenfunctions near 
the boundary. 


1. Introduction and main results 


1.1. Steklov problem. Let Obea compact n-dimensional manifold with C°° boundary dfl 
and unit exterior normal u. We consider the Steklov eigenvalue problem 

A u\(x) = 0, x G fl, 

d„u\(q) = A u x (q), q G <90. (1-1.1) 

The solutions u\ G C ,00 (h2) are called Steklov eigenfunctions. Let 'fan ■ C(O) — > (7(911) de¬ 
note the boundary restriction map. The boundary restrictions of the Steklov eigenfunctions, 
:= jdnU\, are the eigenfunctions of the Dirichlet-to-Neumann (sometimes also referred to 
as the Poincare-Steklov) operator P. It is well-known that P is a first-order homogeneous, 
self-adjoint, elliptic pseudodifferential operator P G T' 1 (cAT), where P — a/— G T°(<90) 
(see hsf Consequently, the Steklov spectrum, which coincides with the spectrum of P, 
consists of an infinite sequence of eigenvalues A j with Xj —> oo as j —> oo. 

Recently, there has been a significant interest in the study of the nodal geometry for the 
Steklov problem, see IBL.i/e. WZ J. However, these results were concerned with the nodal sets 
of the Dirichlet-to-Neumann eigenfunctions tp\, i.e. the boundary nodal sets. The present 
paper focusses on the interior nodal sets, i.e. the nodal sets of the Steklov eigenfunctions U\. 


With the exception of the recent preprint [SWZ] (see Remark 1.2.3), very little has been 
previously known on this subject. We also note that Courant type bounds were earlier 
obtained on the number of the interior (in any dimension) and boundary (in dimension two) 
nodal domains of Steklov eigenfunctions fjKSl lAMl fKKPj. see also EE section 6.1]). 


1.2. Main result. Motivated by a celebrated conjecture of Yau |Y1]IY2] , it has been recently 
suggested in |GPl Open problem 11 (i)] that the ratio between the (n — l)-dimensional 
Hausdorff measure of the nodal set of a Steklov eigenfunction u\ and the corresponding 
Steklov eigenvalue A is bounded above and below by some positive constants depending only 
on the geometry of the manifold. Our main result below gives a positive answer to this 
conjecture for real-analytic Riemannian surfaces. Note that, in this case, the 1-dimensional 
Hausdorff measure of a nodal set Z Ux is simply equal to its total length, which we will denote 
by C{Z UX ). 
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Theorem 1 . 2 . 1 . Let LI be a real-analytic compact Riemannian surface with boundary. Then 
there exist constants Cj = Cj(Q) > 0, j = 1 , 2 , such that, for any Steklov eigenfunction u\ 
with an eigenvalue A > 0, the total length of its nodal set Z Ux satisfies: 

CiA < L{Z Ux ) < C 2 A. 


This result may be viewed as an analogue of the Donnelly-Fefferman bound for the size of 
the nodal set of Laplace eigenfunctions on real-analytic manifolds B3E1- 

Remark 1.2.2. To illustrate the nodal bounds in Theorem |1.2.1 consider the Steklov eigen¬ 
functions on a unit disk corresponding to the double eigenvalue A 271-1 = A 2 n — n, n — 

1, 2,.They can be represented in polar coordinates by u{r6 ) = r n sin(n0 + a) for some 

a G [0,7r/2] . The nodal set of each such eigenfunction is a union of n diameters, of total 
length 2 n. 


Remark 1.2.3. Our methods here are specific to the case of real-analytic surfaces. In higher 
dimensions, Steklov eigenfunctions on 0 are much more complicated and their nodal struc¬ 
tures are not well-understood. In the case of a general smooth n -dimensional manifold with 
smooth boundary, a non-sharp lower bound 

i i „ 2 —n 

\Z\\ > c\^r 


for the volume \Z\\ of the interior nodal sets has recently been proven in [ SWZj . It can be 
viewed as the Steklov analogue of the lower bounds on the size of the nodal sets of Laplace 
eigenfunctions obtained in [ GMl fSZ . HS . M l 1HW] , 

After the first version of the present paper was posted on the archive, an upper bound of 
order A 3 / 2 on the size of the nodal set of Steklov eigenfunctions on surfaces was obtained in 
[ Zh2 ] using a quite different approach. While this bound is not sharp, it is valid for arbitrary 
compact smooth surfaces with boundary. More recently, while the present paper was under 
review, the sharp upper bound in Theorem 1.2.1 has been extended in Zli3l to real-analytic 
Riemannian manifolds of arbitrary dimension. 


1.3. Sketch of the proof. For any closed manifold with boundary Steklov eigenfunctions 
U\ decay rapidly into the interior of fl |HLl GPPSj . In order to analyze their nodal lengths, 
we must therefore consider a neighborhood of the boundary and its complement separately. 
The idea is to use quasimodes near the boundary and frequency function techniques in the 
interior. 

We begin, in section 2, by constructing quasimodes for the Dirichlet-to-Neumann operator 
P. In Proposition 2.3.4 we show that trigonometric functions are approximate eigenfunctions 
for P, up to an exponentially decaying error term. The proof uses the C u surface assumption 
in a very strong way: from |GPPS] . we conclude that the error term decays faster than any 
polynomial in A, and the assumption of analyticity allows us to improve the decay of the 
error term to exponential. From Proposition 2.3.4| and some fairly standard linear algebra 
techniques, we show in Lemma |3.1.1 that the boundary Steklov eigenfunctions ip\ can be 
approximated by trigonometric functions ifx plus an error f\ which decays exponentially in 
A in any C k norm. 

We would like to use the quasimode approximation to estimate nodal length of a Steklov 
eigenfunction u\ near the boundary, but from the example of the annulus (see Example 


1.4), we know that it is possible for u\ to itself be exponentially small near one or more 



































NODAL LENGTH OF STEKLOV EIGENFUNCTIONS 


3 



FIGURE 1. Illustration of the proof. The domain D is the annulus between 
the boundary components dD 1 and dD 2: denoted by solid lines. The dotted 
circles mark a neighborhood of dD x in D and a neighborhood of dD 2 in D c . 
The two sets U C U are as illustrated. 


boundary components. Indeed, this seems to be the generic situation. Near these “residual” 
boundary components, the error in the quasimode approximation can be larger than the 
quasimode itself and hence the approximation is not particularly useful. Our proof avoids 
this difficulty via harmonic extension of the Steklov eigenfunctions across residual bo undar y 
components. By controlling the C ,A -norm of the extended eigenfunctions (see Lemma 4.3.1), 
one can effectively treat a neighborhood of the “residual” boundary components in the same 
way as the interior of D. 

To illustrate the proof, consider a Steklov eigenfunction u\ on the annular domain D in 
Figure[l] Suppose that dD x is a non-residual, i.e. “dominant”, boundary component for u\ 
and that dD 2 is a residual boundary component. In section 3.2, we extend our boundary 
quasimode approximation given by Lemma 3.1.1 into the interior. The exponential decay of 
the error means that the approximation is effective in a A-independent neighborhood of dD x . 
A direct comparison of nodal sets (section 4.2) then gives us upper and lower bounds on the 
nodal length of U\ in a small neighborhood of dD x , denoted by the dotted line in Figure [TJ 
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To treat the interior and a neighborhood of the residual boundary components, we extend 
u\ (section 2.4) to a small A-independent neighborhood of dD 2 , with boundary denoted by 
the dotted circle. In this neighborhood, U\ is still exponentially small in A (section 4.3). Now 
let U C U be two sets as in Figure [l] In section 4.4, we use standard frequency function 
techniques to bound the nodal length of u\ in U from above by the “renormalized Almgren 
frequency function” of the larger domain U. Then, using the quasimode approximation near 
the boundary and the fact that a portion of dU coincides with dDi, we bound this frequency 
function by a multiple of A. A covering argument extends these upper bounds to all of D 
outside a neighborhood of dD 1 , including a neighborhood of dD 2 . Combining these upper 
bounds with the two-sided bounds near dD\ completes the proof. 


Remark 1.3.1. A particularly novel aspect of the problem is the exponential decay into the 
interior of the eigenfunctions, namely that for some positive constants r and C depending 
only on the geometry of f2, 

\u x (x)\<Ce~ Td ^ dn)x . (1.3.2) 


This is an immediate consequence of Lemma 3.2.9 and the exponential decay of the interior 
quasimode approximations. A similar estimate holds for the derivatives of order k > 1 
of u\, provided we multiply the right-hand side by X k . Our problem therefore resembles 
the question of estimating the size of nodal sets in forbidden regions for eigenfunctions of 
Schrodinger operators (see |HZZ1 CT|). A different method of proving (1.3.2) has been 
communicated to us by M. Taylor [Ta2] . 


Remark 1.3.3. It also follows from our results 


surfaces, the error term in the eigenvalue approximation of jGPPS 


see Proposition 3.1.3|) that for real-analytic 

Theorem 1.4] in fact 


decays exponentially as the index of the eigenvalue increases. In particular, for a simply 
connected real-analytic surface 0, there exists a constant r > 0 depending on the geometry 
of Q such that 

2nj 


- X 2j+1 


+ 0(e 


-i-j) _ 


Length(<9G) 


+ 0(e~ TJ ). 


This improves the error estimate 0(j 00 ) obtained independently by Rozenblyum and Guillc- 
min-Mclrose (see |Rol lEd]l for smooth simply connected surfaces. 

1.4. Example: Steklov problem on an annulus. To illustrate several features of Steklov 


eigenvalues and eigenfunctions used in the proof of Theorem 1.2.1 outlined in section L3 
consider an annulus A{ l,e) with inner radius e and outer radius 1. 

First we compare its Steklov spectrum to the spectrum of a union of two disks, of radii 
1 and £, which consists of the double eigenvalue zero, as well as the double eigenvalues k 
and k/e for each k e N. Let us show that the difference between these eigenvalues and 
the corresponding eigenvalues of A{l,e) is in fact exponentially vanishing in k. Indeed, as 
computed in [D] (see also |GPl Section 4.2]), the nonzero Steklov eigenvalues of the annulus 
are given by the roots of the quadratic polynomial 


p k (cr) = a 2 - ak 


£ T 1 \ (1 T £ 


-2k 


1 — £ 


2k 


+ -k 2 , 
£ 
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for each k G N (each root corresponds to a double eigenvalue). We could compute the 
eigenvalues directly, but it is easier to compare Pk(cr) to the polynomial 


Pfc(cr) = cr 2 - ak 


£ + 1 


+h 2 , 

£ 


which has roots k and k/e. In fact Pk(cr) — f(k,e)a = pk(cr), where 

'e + l\ 2e 2k 


f(k,e) = k 


1 — £ 


2k ' 


Since f(k ,£) = 0(k£ 2k ), a straightforward calculation shows that the roots of Pk(cr) differ 
from k and k/£ by 0(k £ 2k ) as well. The error is exponentially decreasing in k, which agrees 


with the eigenvalue approximation result in Proposition 3.1.3 


Steklov eigenfunctions on surfaces also exhibit exponential decay at certain boundary com¬ 
ponents, which motivates our definition of “dominant” and “residual” boundary components 
(see Definition 4.1.1). To illustrate this, consider an eigenfunction u a . of the annulus .4(1, e) 


corresponding to an eigenvalue <jj = k + 0{k£ 2k ). As computed in |GP] , 

k 


U r . 


,(r,e) = c k 


r k + 


a 


J „—k 


k + a,. 


T(kO), 


where T(k6 ) is a unit norm linear combination of cos (kd) and sin(k9) and Ck is an appropriate 
normalizing constant. Since k — Oj = 0(k£ 2k ), we see that Ck ~ 1, that the L 2 -norm of u <Jj 
on {r = 1} is roughly 1, and that the L 2 -norm on {r = £} is 0(£ k ) - which is exponentially 
decaying in k. So, in this case, the circle of radius 1 is a dominant boundary component, 
while the circle of radius £ is a residual boundary component. A similar analysis holds for 
the eigenfunctions corresponding to eigenvalues of the form k/£ + 0(k£ 2k ), with the roles of 
the boundary components reversed. 

Acknowledgments. The authors are grateful to M. Sodin and M. Taylor for useful discus¬ 
sions, and to the anonymous referee for many helpful suggestions on the presentation of the 
paper. The research of I.P. was partially supported by NSERC, FRQNT and Canada Re¬ 
search Chairs Program. The research of D.S was partially supported by NSF EMSW21-RTG 
1045119. The research of J.T. was partially supported by NSERC and FRQNT. I.P. and 
J.T. were also supported by the French National Research Agency project Gerasic-ANR-13- 
BS01-0007-0. 

2. Steklov problem for real-analytic planar domains 


2.1. Koebe uniformization. We first prove Theorem 1.2.1 in the case where Q is a real- 
analytic planar domain and explain the fairly minor modifications needed to treat the general 
case of compact Riemann surfaces with boundary in section 4.5 The two-dimensional case 


is quite special. Indeed, unlike the higher-dimensional case, for surfaces the Steklov operator 
P : C°° (<9f2) —y C°° (dfl) agrees with the square root of the Laplacian on the various boundary 
components modulo a smoothing operator. Upon reparametrization by arclength, the latter 
operator is consequently (modulo smoothing) just the Fourier multiplier M acting on the 
component boundary circles. Although not necessary, one can see the former directly by 
applying conformal mapping. Indeed, by the Koebe uniformization theorem |HSj . one can 
conformally map 9 to a planar domain D, where D is the disk of radius 1 with a finite 
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number of interior disks removed. In the cases where is simply-connected, this reduces to 
the Riemann mapping theorem. When dfl is real-analytic, it follows by C u regularity up 
to the boundary in the associated Dirichlet problem for the Green’s function [ MN ] that the 
conformal map / : D —>■ extends to a real-analytic map of D to fh Moreover, / maps 
boundary to boundary in a univalent fashion. In particular, the induced boundary restriction 
fgci : dD —» dD is a (W-diffeomorp hisrri. 

The boundary dD consists of a union of circles, which we denote by dD i, dD 2 ,..., dDk, 
with radii pi,...,pu and centers Ci,..., c*, respectively. The corresponding boundary com¬ 
ponents of are denoted by dDj := fdn(dDj)-, j = 1, 2,..., k. We let 0 3 G [0, 27r] be the usual 
angle coordinate on dDj for each j, and let the arc length coordinate q 3 be pj9j. Let q be 
an arc length coordinate on dD which coincides with q 3 when we restrict attention to dDj. 
Finally, let dq and dq 3 be the measures on dD and dDj r espect ively induced by the Euclidean 
measure dx on M 2 . As in Pi. the Steklov problem in ( |l.l.l[ ) is conformally mapped to the 
problem 


Av\(x) = 0, x G D, 

d v V\{q) = A gv\(q), q G dD, (2.1.1) 


with g = df(q) = \ f'(q)\ analytic for q G dD and g ^ 0. By conformal mapping, v\ = f*U\ - 
so without loss of generality, it suffices to work with the conformal model (2.1.1) and we will 
do so here. We will abuse notation somewhat and denote the corresponding Steklov operator 
in the conformal model (2.1.1) by P : C°°(dD) —>• C°°(dD ) as well, and its eigenfunctions 
by ip\. We henceforth make the normalization that 


||</k\||L 2 (aD) - 1- 


( 2 . 1 . 2 ) 


Since we consistently work with the model case in ( |2 .1.1 ), this should not lead to confusion. 
Given r = (r 1; ..., t*,) G (0, l) fc , we complexify dD to a Grauert tube (ie. an annulus) 


dD f ; = {(Bj = 1 e i0j ~ Pj ; ( Pj ,9j) G [0,t,) x [0,2tt]}. 


We choose r 3 > 0 here to be the analytic modulus of g\on :i G C“{dDj)\ that is, the maximal 
tube radius for which g has a holomorphic extension to ( dDj )^. 

We also note that the length Lj of dflj = fgn^dDj) is J 9 D .g(qj ) dqj, and define a new 

coordinate on dDj by Sj(qj) = / 0 ?J g(q'j) dql G [0, Lj]. Let s be the corresponding coordi¬ 
nate on all of dD. Note that since g is analytic and strictly positive, this is an analytic 
reparametrization with analytic inverse. 


2.2. Potential layer formulas and the Steklov operator. We briefly review the char¬ 
acterization of the Steklov operator P : C°°(dD) —> C°°(dD) in terms of potential layer 
operators. This material here is well-known and further details can be found in [ED 
Sec. 7.1]. Here, we assume that C is a bounded domain with C , °°-boundary. Let 
G(x,x') G X?'(R” x M n ) be the ambient free Green’s function for —A = — Ej=i i n 
Consider the single and double layer operators SI : C°°(dD) — > C°°(Q U (M n \ f2)) and 
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VI : C°°(dQ) ^C°°(Q U (M n \ Q)) given by 


Slf(x) = - G(x, q) f (q)dS(q), 

Jan 

Vlf(x) = - [ d v ( q) G{x,q)f{q)dS{q), ieflU(r\0). 

JdCl 


( 2 . 2 . 1 ) 


Corresponding to SI and VI are the boundary layer operators S and N : C°°(dQ) —> C°°(dQ) 
given by 


Sf(q) = 


>dn 


G(q, q')f{q')dS(q l ), 


Nf(q) = —2 / d uW) G(q,q')f(q')dS(q'), q G dQ. 


( 2 . 2 . 2 ) 


I an 


These operators are classical pseudodifferential with S,N G \l/ _1 (<9n), and S is elliptic. 
Given a function v G C°(M n \ dQ) and q G dQ, let v+(q) and V-(q) denote the limits of v(x) 
as x —> q from x G and x G M n \ respectively. The layer potential operators in (2.2.1) 
and the induced boundary operators in ( |2.2.2 ) are linked via the boundary jumps equations 

Slf+{q)=Slf4q) = Sf(q), 


VlU(q) = ±\m + \N}(q), 


q G dQ. (2.2.3) 

Now let P be the Dirichlet-to-Neumann operator for Consider the Dirichlet problem 


A u = 0 on u |an = /. 

A straightforward application of Green’s formula gives 

u(x) = VIf{x) — SI Pf(x ), x G fh 


Taking limits in (2.2.4) from within Q and applying (2.2.3) gives 


sp = ~p 


IV). 


(2.2.4) 


(2.2.5) 


Since S G T 1 (9fl) and is elliptic, it follows from (2.2.5) by a parametrix construction in 


the standard pseudodifferential calculus that P G T^cAl) and is also elliptic. 

In the case where dfl is real-analytic, P is also analytic pseudodifferential in the sense of 
EE] (see also m Ch. 5]), and we write P G 4d(<9f2). To see this, we first note integration 
against G(x,x') in M n is a pseudodifferential operator X with full symbol |£| -2 in the usual 
coordinates. As the subprincipal symbols are all zero, the symbol satisfies Cauchy estimates 
and so, X G v h a ( 2 (M n ). Since dQ is real-analytic, the Fermi coordinates (x,x n ) near dQ are 
also real-analytic. In terms of these local coordinates, ui(x, x n , £, f n ) is thus an analytic 
symbol; hence so is d Xn ax{x,x n ,fjf n ). It is easy to check that the symbols of S and N are 
given, respectively, by J R (cr x | Xn=0 ) df n and 2 f R (( d Xn tr x + if n cr x ) L=o) d£ n respectively, and 
hence are also analytic. So S, N G T“ 1 (cAl). By constructing a parametrix for S in the 
analytic pseudodifferential calculus (see na Ch. 5] for details), then multiplying ( |2.2.5[ ) by 
this parametrix on the left and rearranging, it follows that 

P e 
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In the following, we say that an operator R E \I/*(<9n) is analytic smoothing if its Schwartz 
kernel E C u (dVL x dD), and we write R E ^f°°(dD). 


Remark 2.2.6. In what follows, we use the following notation: given a set X and two non¬ 
negative functions / : X —>• R and g : X —> M, throughout the paper we write / ;$ g if there 
exists a constant C > 0 such that f(x) < Cg(x ) for all x E X. In addition, f & g will mean 
that both / ;$ g and g (S f are satisfied. 

2.3. Quasimodes for the Steklov operator. It is well-known f GPPS . see also |EdlIHLj ) 
that for smooth bounded planar domains, 


P = M + R, Re V-°°(dD), 


(2.3.1) 


where M = ©LMj, with Mj : C°°{dDj) —» C°°(dDj) the Fourier multiplier defined by 


M. e ^imsj/ L j _ 


27 r|m| 

~T~ 


^imsRLj m e 


Note that M = \Z—Ag D . As in [ GPPSJ . let A(a) be the sequence {0, a, a, 2a, 2a ,... } and 
A(a i,..., CKfc) be the union of these sequences, with multiplicity, arranged in the appropriate 
order. We write 




\Li 


Lk ) 


(2.3.2) 


which is the nth eigenvalue of the operator M. As a consequence of (2.3.1) and [GPPS], the 
eigenvalues of P consist of the sequence 

A n = /U n + 0{n 00 ), n E 

It follows that X n ~ n as n —> oo. We abuse notation somewhat and denote the corresponding 
orthonormal basis of eigenfunctions of P on dfl by so that each is the restriction 

of u\ n to dQ. Although the operator M acting on dD is non-local, it clearly maps each of 
the boundary components dDj to itself; that is, for any / E C°°(dD ) with supp / C dDj, 

Mf = Mjf. 

Denote the L 2 -normalized eigenfunctions of M on the boundary component dDj by e m j(sj), 
where 


e m As j ) = e 2 * ims i' L i, Sj E [0, Lj). 


(2.3.3) 


We also let e m j(s) be the function which is e m j(sj ) on dDj and 0 on the other boundary 
components, with an analogous definition for tp\ n j- 

The following lemma is central to the proof of Theorem 1.2.1 and shows that when fl is 
analytic, the functions e m j are quasimodes for P to exponential error in m. 

Proposition 2.3.4. Suppose that D is analytic, with D, the conformal map f : D Q, 
and the Dirichlet-to-Neumann operator P : C°°(dD) —y C°°(dD) as in (2.1.1). Then: 

(1) The remainder R : C°°(dD) —» C°°(dD) in (2.3.1) has Schwartz kernel R E C u (dDx 
dD) with analytic modulus at least t > 0 in each variable separately. Here, t depends 
on the analytic modulus of the conformal multiplier g E C u and the geometry of D. 
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(2) As \m\ -A oo ; 

II Pe m j(a) - ^AAe m j( S )\\ CHgDj = O k ( e - 2 ’ f H/%). 

In the last line, O k (e~ 27TT ^ m d L J) denotes a non-negative term that is < C k e^ 27r d m \/ L 3 
with constant C k > 0 depending on k. 

Proof. Our proof uses the boundary integral equation for ipx- For x , x' G M 2 , we continue to 
let G(x, x ') be the free Green’s kernel for — A K 2 = — d 2 — d |2 ,, which is — log \x — x'\. Using 


(2.2.5) (see also 


oj), the functions \ satisfy the boundary jumps equation: 


A</W(g) = A n / G(q,q')ip Xn (q') 9W)dq' ~ / K(q, q')<p\(<f) dq ’, (2.3.5) 


'dD 


with 


q ) &v{x')G{x, X ) | ( x,x')={q,q') 


JdD 

1 d^ x ’)\x - x‘ 
47r \x — x‘ 


./12 


/12 


\(x,x’)=(q,q')- 


Change variables in the first integral to integrate in terms of s': 


aVA„(<?(s)) = A n / G(q(s), q(s , ))(p\ ri (q(s')) ds' — / K(q(s), q(s'))(px(q(s’)) dq(s'). (2.3.6) 




I dD 


We now claim: 


Claim 2.3.7. 


2^ a „ = \ n G<px n ~ KKiP\ n ~ K 2 (p Xn , 


where AT, K 2 G v l/ a 00 (9Zi)) and G = ®j =1 Gj with Gj : C°°(dDj ) —>• C°°(dDj) given by 


h{s')ds'. 

Recall that is the space of operators with C u Schwartz kernels. We also claim: 

Claim 2.3.8. If we extend Mf 1 to act on C°°(dDj ) by letting it be zero when acting on 
constants, 


( Gjh)(s) = / log \e 


2ivis/Lj _ ^ivis'/Lj | 


Mr 1 = 2Gj. 


Assuming both claims, we now prove Proposition 2.3.4 


The claims combine to show 


<p Xn = A n M Va„ - 2\ n K 1 (p Xn - 2K 2 <px n = M 1 Pip Xn ~ 2\ n K 1 p Xn ~ 2 K 2 <p Xn . 

Given / G C°°{dD ), we can write / = Y%=i a kV\ k with a k = 0(fc -00 ). By linearity, the 
boundedness of K\ and K 2 , and the rapid decay of a k and of a k X k , it follows that 

M~ l Pf -f = 2K\Pf + 2 K 2 f. (2.3.9) 


Since K \, K 2 G T a °°(dD ) and P G \l/^( dD ), standard theory of analytic pseudodifferential 
operators m Ch. 5] implies that K = 2K\P + 2K 2 G Af a °°(dD), with 

M~ l Pf — f = Kf , or equivalently (P — M)f = MKf. 


(2.3.10) 
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Now, we choose f(s) = e m j(s). For any r < Tj, by contour deformation we have 
Ke md (s)= [ K(s, Sj )e 2 ^^d Sj 

i 

= ( [ K(s , Sj ± if)e 2nimSj/Lj dsj'\ = 0( e - 2 *M*7^-) > (2.3.11) 

' 'dDi ' 


< dD, 


where the choice of ± is determined by the sign of m. The same is true for each (d k Ke m )(s), 
so for all k, 

\\ke mJ \\ ck = O k (e- 2 ^ L >). 

The same is therefore true for ||/Ce m j||#fc for any k. Since M G '$! 1 (dD) is a (classical) 
pseudodifferential operator, it is bounded from H k —y H k for each k. From this and 
Sobolev embedding we conclude that \\MKe m j\\ck = Ok(e~ 2n l m l r °/ L i) for any k. Combining 
this with (2.3.10) completes the proof of the proposition. 


□ 


We now prove both claims. 


Proof of Claim 2.3. 7 . Assume that q(s) G dDj. The portions of the integrals in (2.3.6) over 
dD \ dDj can be absorbed into K\ and Jfo, since d(dDi, dDj) > 0 for i ^ j and therefore 
the integrands are real-analytic functions of s'. Using the symbol = to denote equivalence 
up to terms of the form (X n Ki + K 2 )(f\ n with K\ and K 2 G we have 


2 < / ? A„(fo( s j)) — X % 


G{qj(sj), qjis'j^pxAqjis'j)) ds'j 


ldD.j 


K (0j(sj), QjWj^Vxiqjis'j)) dqj(s'j). (2.3.12) 


For notational simplicity, we now suppress all j subscripts. The last term in (2.3.12) can 
be absorbed, because the normal derivative K(q(s),q(s')) is itself real-analytic in (s, s') G 
[0,L] x [0,L]. To see this, observe that analyticity away from the diagonal s = s' is clear. 
To establish the analyticity near s = s' we insert the Taylor expansion q(s) = q(s') + (s — 
s')d s >q(s') + ^^1 q(s) — q(s')\ 2 v(s') + [ar(s, s')d s iq(s') + aN(s, s')v(s')} into the formula for 
K(q(s), q(s')) in (2.3.5). Here k(s) denotes the curvature at the point s. The result is that 

a N (s,s') 


rm ,«w> = -A(^ + 


?( s 0l" 


and so, (s, s') h- > K(q(s), q(s')) G C^QO, L\ x [0, L\) since ajy is real-analytic with a^(s, s') = 
0(|s — s'| 3 ) as s -> s'. 

Finally, all we need to do to complete the proof is to show that 


1 

27T 


log I q(s) - q(s')\ 


— log \e 2nis / L — e 2 ' iris '/ L 
27r 


1 lor 

27T ^ |g27rz5/L _ e 2nis'/L\ 


is analytic. However, the reparametrization q —y s(q) is analytic with analytic inverse, 
and Taylor expansion of the quotient on the right-hand side shows that it is analytic and 
nonzero at s = s' (cf. ( |TZl formula (8.8)])). Therefore, its logarithm is also analytic there. 
Analyticity away from s = s' is automatic. □ 
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Proof of Claim 2.3.8 . By scaling, it suffices to prove the claim for j — 1 and L\ = 2n. We 
will show 


1 
7r 


-2?r 


log \e ls — e 


ls '' e ims ds' = -~e ims , meZ\{0}. 
\m\ 


(2.3.13) 


The claim follows since {e irns } ne z is a Fourier basis. To prove (2.3.13), consider the harmonic 


function on the disk whose Neumann data is e ims . On the one hand, this is 

’ \m\ 

other hand, the Poisson kernel for the Neumann problem is — - log \x — x'\, so 


On the 


,|m| 


-2it 


gtms _ 


\m\ 


-log I re ls — e 

7r 


is | gims 


ds'. 


Restricting to r = 1 shows (2.3.13). 


□ 


2.4. Harmonic extension of interior eigenfunctions across the boundary. It is well 

known that a harmonic function on a domain can be extended harmonically across a real 
analytic bound ary (see [EM])- The following lemma, which will be used in the proof of 
Theorem 1.2.1 provides an explicit estimate on the (W-norms of the harmonic extension of 
a Steklov eigenfunction u\ across a boundary component Dj in terms of the L 2 -norm of its 
trace Lp\ on dDj. 

Proposition 2.4.1. Let u\ G C U (D) be a Steklov eigenfunction. Then, for every connected 


component dDj of the boundary, and 0 < e < r with r defined in (3.2.1), there is a harmonic 
continuation ofu\ (also denoted u\) across the boundary component dDj to an annulus A~j{e) 
of width e. Moreover, for k = 0,1, 2,... one has the exterior estimate 


ll' U A||cfe(A+(e)) < C kjE \ k+1 e X£ \\ip x \\ l 2 (8D 3 


Proof. From Proposition 2.3.4, we know P = M + R where R G 'i , a °°(dD). We use this 


characterization of P to estimate Fourier coefficients of <p>\ \go :i along the boundary circle 
dDj. To simplify the writing in the following we assume that Lj = 2ti and that the center 
of the disc Dj is the origin. Also, we abuse notation somewhat and simply write <^a(s) f° r 


<Px(qj(s)) in ( |2A2| -( [2A4| ) below. 

Then, for k G Z \ {0}, 

•fix(k) = {e isk , ip\) = T= T<e“, (P - R)<p,) 

= w&'' k ^ }+ isi 0(e ^ H)||w|1 = wr (k)+ isi o(e ^ H)|1 ^ 11 - 


From (2.4.2), we observe 


(2.4.2) 


(2.4.3) 


(\k\ - X)^(k) = 0(e- T l‘l), 

and hence, 

fi x {k) = 0(e~ T \ k I), |fc|>A + l. (2.4.4) 

It follows that ip\ holomorphically continues to the strip [0,27 t] + i[—£,e\ and so ip\\ qd 6 
holomorphically continues to an annular neighbourhood of dDj of any width £ < r. In terms 
of the parametrizing coordinates s G [0, 2n\, one holomorphically continues (p\ to the annulus 

Aj(e) :={igC;i = e lw , w G [0, 27r] + i[—£, e]}. 
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Without loss of generality, we assume here that the set where Imw < 0 (denoted by Aj (e)) 
is the part of the annulus lying inside the domain D and Irn w > 0 (denoted b y Af( e)) is 
the part lying outside. Note that £ > 0 is independent of A. It also follows from (2.4.3) that 
the holomorphic continuation, of the boundary Steklov eigenfunction p x satishes 


1 ^ 0*01 < C e e 


A|Im w\ 


\WA\l? {dDj), 


x = e 


(2.4.5) 


with an appropriate constant C £ > 0. 

Now we need to construct the harmonic continuation of the interior Steklov eigenfunction 
u\ E C W (D). By Green’s formula, 


u x (x) = (Vlp x )(x) - X(Slgp x )(x), x E D. 


(2.4.6) 


Since the harmonic continuation of both terms in (2.4.6) is carried out in the same way, we 
consider here the second term and then just indicate the minor changes necessary to deal 
with the first one. Since the multiplicative factor A is irrelevant to the continuation, we just 
consider 


(Slgip x )(x) = — 


log \x - q | g(q) <p x (q) dq + Rip x (x), x E A - (5) 


(2.4.7) 


’on. 


where Rp x = A. f an , an log \x — q\ g(q) p x (q)dq. That Rp x (q) harmonically continues to 

x E is immediate since the kernels G(x,q ) with (x,q) E Aj(s) x dD^,j ^ k are 

supported off-diagonal and so G(x, q) is harmonic in each variable. Thus, it is enough to 
analyze the first term on the RHS of (2.4.7). 

The arclength parametrization of dDj is given by q(s) = e ls E dDj, s E [0, 27r]__We let 
©a (s) := /o g(q(s'))(p x (q(s'))dq(s'). First, by application of Green’s formula in (2.1.1) with 
harmonic measure 5qd p it follows that g(q(s'))(f X (q(s'))dq(s') = 0 and so, in particular, 
©a G C^QO, 27 t]). Moreover, since 


\9x( k )\~\ k \ \9<P\(k)\, 


it has the same holomorphic continuation properties as gip x (s) (with continuation to the 
strip [0,27 t] +*[—£,£]). Then, making the definition 

$ A (e“) := 0 A (s), sG[0,2tt], 


it follows that 4 >a E C u (dDj ) holomorphically continues to the £-annulus Aj(e) and satishes 
(2.4.5) there. Then the RHS of (2.4.7) can be written in the form 


2n 

log \x - e ls \ d s ($ x (x) - $ x (e ls )) ds + Rp x {x). 

Here, we continue to view x = e lw E C as a complex variable. Then, writing 2 log x — e ls = 
log (a; — e ls ) + log(x — e ~ ls ), an integration by parts with respect to d s gives 

1 r 2n — * Ap is ) 

Sl(gp x ){x) = —Re / —--- -ie ls ds + Rp x (x); x E Aj(e). (2.4.8) 

Z7T j o x — e 

By Taylor expansion and Morera’s theorem, for each s E [0, 27r], the function f X)S {x) : = 

^ ^ ias same analyticity properties as p x (x) and consequently, fx, s extends to a 
holomorphic function for all x E Aj(s). Thus, the first term on the RHS extends harmonically 
to x E Aj(e). 
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Also, we note that for each x G Aj(e) the function Fx, x (s) = 4,x ^_gi^ e ^ extends holomor- 
phically to s G [0,27r] + By Cauchy’s formula, we can therefore deform the contour 

of integration into Aj(e) by letting s H> s — iSo for any 0 < 5 0 < £■ Undoing the integration 
by parts in (2.4.8) then implies that we can write the harmonic continuation formula for the 
single layer term in the form: 


Sl{g<p\)(x) = Re 


-2 IT 


log \x—e ls+So \ g(s—ido) (px(e ls+ °)ds+R(p\(x), xeA^(e), (2.4.9) 


where log in (2.4.9) d enotes the principal logarithm (we note that \x — e ls e s °\ > 0 when x G 
Aj~(e)). The formula (2.4.9) is useful for estimating the harmonic continuation of Sl(gip) into 
the exterior of D (where it can blow-up exponentially in A). In deed, since j log(x — e ls+s °)| = 
0(1) uniformly for (x, s) G A^(S 0 /2) x [0,27r], it follows from (2.4.9) and (2.4.5) that 


\Sl(g(fx)(x)\ < e XSo \\ip x \\L*(d Dj ), 


(2.4.10) 


uniformly for x G A+W2). 

The double layer term T>l((p\)(x) is analyzed similarity. In this case, we again rewrite the 
integral in complex form. Given the parametrization [0, 27 t] 9 s t— y e ls = q(s) of dDj, the 
unit normal is u(q(s)) = e ls and so, 


dv(q(s))G(x, q(s )) = d„( q ( 9 )) log \x - q(s ) | = Re 
As a result, we can write in complex form: 

q(s. 


x — e l 


%»)(*) =Re [ ■ ( -^ ) ds. 


(2.4.11) 


Next, writing ^Ayc) s log(x — q(s)) = one integrates by parts in s with the effect of 

replacing the 0(\x — g(s)| _1 ) singularity in (2.4.11) with a log \x — q , (s)|-type singularity. In 
the process, one creates an extra factor of A coming from differentiation of (p\{q{s)). To see 
this, we note that from Proposition 2.3.4, P = ®jMj + R with R G T” 00 and Ptp\ = \<p\. 
Since Mj(e l7nSj ) = m \ e mS:i . for all m, it follows from Fourier expansion that 

WdstpxWtfidDj) = \\Mj<px\\ L 2( dDj ) = ||(R - R)<P\\\L2(dD j ) = 0(A)\\(p\\\ L 2( dDj ). 

So, modulo replacing (px(q(s)) with d s tp\(q(s)), the analysis proceeds as in the single-layer 
case and \T>l(tp\)(x)\ also satisfies the bound 


\Vl(ip x )(x)\ < Ae A5o ||^ A || L 2 (5 ^.), 


(2.4.12) 


uniformly for x = e lw e A+(5 0 / 2). 

Finally, we note that for k > 1, \\Sl(gipx)\\c k an d ||X , /( < ^ A )||c 7 fc are estimated in the same 
way, except that one must first integrate by parts k times in s to compensate for the k - 
derivatives in x. This creates an extra multiplicative factor of X k . As a consequence of (2.4.6)- 


(2.4.12), we have proved Proposition 2.4.1 


□ 
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3. Approximation of Steklov Eigenfunctions 

In this section, we construct approximations to our Steklov eigenfunctions, which we will 
use in our study of their nodal sets. First we use Proposition 2.3.4 and our functions e m j to 
construct these approximations on the boundary. Then we show how to construct quasimodes 
for Steklov eigenfunctions in the interior. 


3.1. Approximation of boundary eigenfunctions. As a first step, Proposition 2.3.4 


allows us to show that the boundary Steklov eigenfunctions <px n can be approximated up to 
exponentially decaying error by linear combinations of the quasimodes e mj -. More specifically, 
one has the following 

Lemma 3 . 1 . 1 . There exist functions ifx„ and f\ n in C°°(dD), n = 0,1,2,..., such that 
ipx n = 0 x n + 1'k f or eac/i n ! an d furthermore: 

(1) Let Lpx ri ,j ■— 0aJ 8Dj, and suppose n is sufficiently large. Then for each boundary 
component dDj, either ifx n ,j = 0 or there exist m n j e Z as well asb n j j+ andb n j_ e M. 
such that 

J, . _ U ImmnjSjILj , > -2irim nij s d /Lj 

rA n ,j — u n,j,+e , 

with 

I \ 27 Tm n,j | , 271 J rj T , 

I An- j—^-\ < —, L maxjLi,..., L k \. 

(2) For each i € Z + , there exist constants C^> 0 and r 0 > 0 such that 

\\f\n\\c*(dD) < CeX^e T ° Xn . 

Remark 3.1.2. Since <px n are L 2 -normalized, for sufficiently large n, there exist boundary 
components for which ifx n ,j are n °t too small - in particular, large compared to \fx n \dD y 
These boundary components will be called dominant for ux n (the precise definition follows 
in section 4.1). Those components on which ifx„,j are small, or vanish identically, will be 
called residual. 


Another consequence of Proposition 2.3.4 is eigenvalue asymptotics with exponentially 
decaying error: 

Proposition 3.1.3. There exist C > 0 and iq > 0 such that |A n — p n \ < Ce~ Tin . 


The proofs of Lemma 3.1.1 and Proposition 3.1.3 are relatively standard, but involve some 
technical linear algebra and are thus deferred to the Appendix. 


Remark 3.1.4. From Proposition |3.1.3] for each eigenvalue X n with n sufficiently large, there 
exists at least one boundary component dDj and an integer m n j for which 

3.2. Interior quasimodes. In this section, we extend our approximate boundary eigen¬ 
functions ifx n to the interior and show that the extensions are good approximations to the 
interior Steklov eigenfunctions. Recall that 


0A, = 


k 

E 


(b,, , . ( i-n/ 7 -. _|_ b n j . 
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We will extend each to an almost harmonic function u n j e C°°(D), then sum them to 
create our global interior quasimode u n . 

As a hrst step, we dehne the rate constant 


r : = nun 

j=l,...,k 


2vr 

To, -j-r 
W 


(3.2.1) 


where r 0 is the exponential rate for the approximation in Lemma |3.1.1| and r is the minimal 
analytic modulus in Proposition |2.3.4 


To construct the needed approximately harmonic functions, let Sj : dDj —> [0, Lj\ be the 
arc length function (with metric g) along the boundary circle dDj. Without loss of generality, 
we assume here that the circle dDj is centered at the origin. Identifying dDj with M/27rZ 
via the parametrization 6 : [0, 27t] —y dDj , we let 


r ,- := 


mm s 

0S[O,27t] 


'M 


Since Sj is monotone increasing, each of the constants Tj > 0. We let sj(z) be the holomorphic 
continuation to the strip S T := [0, 27r] +i(— t,t) and put Nj := ma xg + ^ e s T |<9|lm sf(9 + i£)|. 
We now £x a constant S > 0 once and for all: 


5 = 5 ( r ) : = “in { r w/2 - NjT 3 , r/2}, 

7=1 ,-,k 


(3.2.2) 


with r > 0 dehned by (3.2.1). By possibly shrinking r > 0 further, we can (and do) assume 
from now on that 5 > 0. 

In the following, we identify holomorphic extensions of 27T-periodic C^-functions to a strip 

iz — 

over [0, 27 t] with an annular neighbourhood of \z\ = 1 via the conformal map z e->■ e L i , 

[0, 27r] + i(—r, t). For each j, let 

e nji z ) '■= < '-.i z G [0,2tt] + i(— t, r). 

Then, natural harmonic extensions to the strip are given by 

u n,j,+ ( z ) '■= Ree nj(-), u n,j,-(z) = Ime^-(z), 0 G [0,27r] +i(-T,T). 

With z = 9 + i£, we can write 

,2n 

■Ti 


u 


n,3, 


[z = e 




/ x sDz) 

u n ,j- [z) = e D J sill 


cos (—m n jReSj(z)), 
'i 

2t r 


(—m n}j Resf(z)), 


where, by Taylor expansion, Re s^ and Im sj are real-analytic functions with 

Ims?(9,«) = «'(»)« + OK 3 ), 


Re4(0,f) = s(«) + O(f). 

As above, since the arclength function s(9) is strictly increasing, 


(3.2.3) 


(3.2.4) 


s\6) > Tj > 0. (3.2.5) 

To dehne the global quasimode corresponding to U\ n we let Xj G (^“(M 2 ; [0,1]) be standard 
cutoff equal to 1 in annular neighbourhood of dDj of width r/2 and vanishing outside an 
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annulus of width r. Then, with x = e tz , and j G {1, k} corresponding to an outer boundary 
component and with b n j t ± as in Lemma 3.1.1 we define the global interior quasimode to be 

Un,j (t) • X j ( X ) ‘ (bn (^) T b n ,j t — U n ,j,— (<z)) 


= Xj\ x ) 


bn,j,+ 6 


<C/ 


/27T 

cos (—m n jRe s “ (2 
W 


)) 


+ b n ,j,~ e 


— ffrmnjlm s^(z) 


,2vr 


sin (—m n jRes^(z))). (3.2.6) 

% 


We note that in the case of an outer boundary component, the point x = e lz with Irn z = 
£ > 0 will lie in an interior collar neighbourhood of t he bo undary circle. In view of (3.2.4), 
for such points, Ims^ ( 2 ) > 0 and so, the quasimode (3.2.6) decays exponentially in n. 
When dDj with j G {1,..., k} is an inner boundary component, we replace e^j(z) above 


with e„ ( 2 ) := e L i 




and replace u n j.±(z) in (]3.2.6b with Ree^ Az) and Ime^ -(z) 


respectively. In this case, the corresponding quasimodes decay exponentially in an interior 
collar neighbourhood where Im 4M < 0. 

The function u n j G C°°(D) is approximately harmonic in D and agrees with the boundary 
quasimode fi>\ n ,j on dDj. Indeed, since A u n j t ± = 0, by Leibniz’ rule for the Laplacian 

A u n j A (xj • (^n,j,+^ n IT T b n j—U n j- ) 

2V xX.j ' ^ (bn,j,+^n,j,+ T b n ,j,-U n ,j,-) T A Xj (^nj,+H?i,J,+ T b r j -ll n j - ) . 

Since the derivatives of Xj are supported in an annular neighbourhood of dDj where r/2 < 
\x\ < r, it follows in view of (3.2.4) that with <5 = 5(r) as in (3.2.2), 

IIAx^njIlc^D) = 0(X n e SXn ) 

and similarity, 


II A x u n ,j\\c*(D) — e SXn ), l G Zb 

Now for each n, we define the interior quasimode 

k 

U n . ^ ) U r 


“n,j- 


(3.2.7) 


(3.2.8) 


J=i 


As we now prove, the function u n in (3.2.8) is the quasimode approximation to the Steklov 
eigenfunction u Xn in a collar neighbourhood of the boundary dD that we seek. 

Lemma 3.2.9. For any fixed £ G Z + , there exist constants Cj(£ ) > 0; j = 1, 2, such that 
Proof. For each n, let 

fn-=u Xn -u n , (3.2.10) 

so that f n is the error in the interior quasimode approximation. To ||/ n ||chD) we no ^ e that 
by construction, u n = i/j\ n on dD, so J n \dD = f\ n , which by Lemma 3.1.1 and the definition 
of 6 shows that 


WfnldDWctOD) < C k X^~ ToXn < CXe~ 5Xn . 


( 3 . 2 . 11 ) 
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Moreover, both u n and f n are approximately harmonic. Indeed, in view of (3.2.7) and (3.2.8), 
we see 


| fn\ \ (D) — || < ClX^e 


To begin, we claim that 


ll/n||c°(-D) — CX n e <5A " 


(3.2.12) 


(3.2.13) 


Indeed, let $ be the solution of A$ = — 1 on D satisfying the Dirichlet boundary condition 
$|qd = 0. Then, by the maximum principle, $ > 0 in D. We define 


g n :=C 0 e-^ + CiA n e- flA "$ >0 


5 A?? 


By (3.2.11) and (3.2.12), f n < g n on dD with A f n > A g n . Similarity, f n > —g n on dD with 
A(—g n ) > A f n on D. Consequently, again by the maximum principle, 

I fn | < \g n \ = 9n on D, 


which proves (3.2.13). 

To treat the cases where m > 1, we note that by standard elliptic estimates ( [Kr] . Theorem 
5.4.1), for any m G Z + even, 


\\fn\\ H rn (D) A C m ( 11A 2 /n||L 2 (_D) + ll/ra||L 2 (D) + H/nlaol 


H m -?(dD) 


J- 


Substitution of (3.2.13) and (3.2.11) in (3.2.15) gives 

D) A Cm (|| A 2 fn\\L 2 {D) + ^n e ^ 


u 


n\\H m (D) 


+ An 2 e 


0- 


(3.2.14) 


(3.2.15) 


Finallly, using (3.2.12) and the fact that \\A. m t 2 f n \\ L 2 ^ ^ ||A/ n || C m/ 2 -i(re), it follows that 

II fn\\ H m (D) = O m ( (An 2 + An 2 + A n ) e~ 5Xn ) = O m (max{An 2 , An} e~ SXn ). 

The proof then follows by an application of the Sobolev lemma which gives 

\\fn\\c e {D) ~ \\fn\\H m (D)i m > t + 1- 


□ 


4. Proof of Theorem 11.2.11 


We now prove Theorem 1.2.1 first we prove lower bounds on nodal length and then 


proceed to upper bounds. Recall that U\ is the Steklov eigenfunction on Q with eigenvalue 
A, and V\ = f*U\ is the corresponding eigenfunction on D. Throughout, we denote the nodal 
length of a function u on a set A by C(u,A). Note that since the Steklov eigenfunction in 
the conformal model is V\ = f*U\, the ratio of C(u\,A) and £(v\, f*(A)) is bounded from 
above and below by geometric constants independent of A and of A C 0; hence we may work 
with the conformal model D without creating problems for our estimates. To simplify the 
writing, we abuse notation somewhat and henceforth denote both the Steklov eigenfunction 
and its image in the conformal model simply by U\. 
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4.1. Dominant and residual boundary components. Our first step is to group the 
boundary components dDj;j = 1, k, into dominant and residual categories. The idea here 
is that since the basic quasimodes approximate actual Steklov eigenfunctions to 0(e~ rA ) 
in C k norm (by Lemma 3.1.1), eigenfunctions that have L 2 -norm less than e _rA along a 


boundary component have no meaningful quasimode approximations 

A, 


n (see section 


We note that since 
2.3h , the bounds in Lemma 3.1.1 in terms of e _T ° n are comparable to 


ones in terms of e~ TAn by possibly adjusting the constant r > 0. We choose to work here 
with the latter. The simple example of the annulus (see Example 1.4) shows that L 2 -norms 
of Steklov eigenfunctions can indeed be exponentially small with ||^a||l 2 ( 9 Dj) ~ e -AC along 
certain boundary circles with some C > 0. Moreover, it is not clear in general that C > 0 
will be smaller than the exponential rate r > 0 appearing in the quasimode approximation 
in Lemma [3.1.11 


DEFINITION 4.1.1. Recall the definition of 5 from (3.2.2). The boundary component dDj 
is said to be dominant for u\ provided 


\ u x\\l 2 (dDj) > e sx ^ 3 . 


Otherwise, it is said to be residual. 


We refer to Example |1.4| for an illustration of the concepts of dominant and residual bound¬ 
ary components. If a surface has a single boundary component, it is obviously dominant for 
each u\. 

A key observation about our dominant and residual boundary components is the following. 
Recall the approximation in Lemma |3.1.1[ where the approximate boundary eigenfunction 


f>x n is a linear combination of trigonometric polynomials on each boundary component dDj 
with coefficients b n j ±. The following is an immediate consequence of our definition of dom¬ 
inant and residual, combined with the error bounds in Lemma [3.1.1 


. 1 

•) > y 


2<5A^ 
— 3 ~ 


Proposition 4.1.2. If dDj is dominant for u\ n , then 

IWI 2 + l^,-r>e-^+0(e- 4 ^ 

as long as n is sufficiently large. 

Remark 4.1.3. We note that even though L 2 -mass can be exponentially small along bound¬ 
ary components, it has recently been proved by a Carleman argument applied to Steklov 
eigenfunctions on the boundary dfl ( jZhll Theorem 1]), that they satisfy quantitative unique 
continuation in any dimension. More precisely, there is a constant C = C(Q) > 0, such that 

11^11^(01,-) > e“ CA Vj = 1,..., k. 
where dDj] j — 1,..., k denote disjoint subsets of <90. 


(4.1.4) 


Note that as A —> oo, the dominant components of the boundary change depending on A. 
However, there are only finitely many configurations of dominant and residual components, so 
by passing to subsequences we may consider each configuration separately. Therefore, with¬ 
out loss of generality, we assume that the boundary components dDj] j = 1,..., P are domi¬ 
nant and t he rem aining ones with j — P + 1,..., k are residual. Given the L 2 - normal iz at ion 
condition (2.1.2), we may also assume that 11 wa| |sz?i A A (at least one dominant boundary 
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component has this property for sufficiently large A, and there are only finitely many con¬ 
figurations). In the following, we denote by Aj(a) an a-width annular neighbourhood of the 
dominant boundary component dDj. 


4.2. Nodal length near a dominant boundary component. We now obtain optimal 
upper and lower bounds for the nodal length of a Steklov eigenfunction u x near a dominant 
boundary component. 

Lemma 4.2.1. Let dDj be a dominant boundary component and let r > 0 and 5(t) > 0 
be defined as above (see 3.2.2 and 3.2.1). Then, there exists a = a(r) > 0 and geometric 
constants Cj = Cj(dD , D)\j = 1, 2, such that, for any Steklov eigenfunction u\ 

Ci A < C{Z Ux n Aj(a)) <C 2 X. 

Proof . We use the decomposition u Xn =: u n j + f n j, with C fc -bounds on f n j as in Lemma 

Rescale all of these fu nction s by multiplying by C > 0, 

for sufficiently 


3.2.9 (note that f n j = f n in Aj(a)). 


where C is chosen so that b 2 j + + b^j_ = 1. By Proposition 


4.1.2 


U<2e^ 


large n. None of the zero sets change, and by our observations about the size of C, we now 
have 

I \fn,j | \c k = O k (e- sx ^ 2+5x -/ 3 ) = Oh(e~ SXn//6 ). (4.2.2) 

For simplicity, assume that dDj is the outer boundary (a similar argument works for the 
inner boundaries). As in section [R2 we use the coordinates (#,£) G [0,27 t] + i(0,d) in the 
5-strip model of the annulus A(a) with complex variable z = 9 + if. Here, £ > 0 corresponds 
to the interior of D where eigenfunctions decay exponentially in X n . A direct computation 
with the quasimodes of (3.2.6) gives 


-2 


^ n,j X n 


| d e U Ut 


12 > 

\ <, e > 3 


(4.2.3) 


We recall from (3.2.4) that lms^(0,£) = s'-(0)£ + 0(£ 3 ) with Tj > s'j{9 ) > Tj > 0 and so 
lms^(0, £) p» £ for 0 < £ < C'i(r) with Cfir) small. Thus, choosing C'i(r) > 0 sufficiently 
small, one can arrange that 

0 < Im sf(9,£) < 2C 1 (r)f j . 

Then, by possibly shrinking C\{r) > 0 in (|4.2.3) further, one can arrange that for 0 < £ < 

Ci(r), 


u n,j + 


-2 


| dgUn^ 


12 O Sl7l ri) j/b _ S\ n /5 


(4.2.4) 


In the last inequality in (4.2.4), we have use the fact that | m n j — A„| < 27t/L. 

Now, £x a = Cfir) and let E Xn be the subset of Aj(a) = {(0,£);O < £ < a} where 
\h n ,j\ 2 < e- 5 ^/ 4 . By formula (4.2.2) , for A n sufficiently large, Z UXn fl Aj{a) C E Xn . But 
using (4.2.4), for A n sufficiently large, 

\d e u n ,j{9,i)\ > CX n e-^ lras ^° > C'X n e~ sx ^ w , (0,£) e E Xn . (4.2.5) 

Using (4.2.2), eigenfunctions are approximated by quasimodes in C' fc -norm to 0(e~ SXn Z 6 )- 
error and so, it follows from (4.2.5) that for the actual eigenfunctions, 

27T m n j C//i 

|^«a„(0,£)| > CX n e~^ > CdX n e~^A\ (0,£) G E Xn . 
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So on each connected component of E Xn , u Xn is either monotonically increasing or mono- 
tonically decreasing in 6. Since for each fixed £, u n j has exactly one zero in each connected 
component of E Xn , so must u Xn for each hxed £. From (4.2.5), it follows by the analytic im¬ 
plicit function theorem that the nodal set Z Ux D Aj(a) is a union of C u curves, one in each 
connected component of E Xn , which are graphs with dependent variable 9 and independent 
variable £. There are precisely 2/i n zeroes of u n j, so we can write Z Ux D Aj(a) as a union 
of graphs of C u functions /i(£),..., f2^, (0 ■ 

It remains to control the geometry of these graphs. To do this, we claim that there exists 
a constant C > 0 such that 


I%gJ < C'\d () u Xn \ in E Xn 


(4.2.6) 


Given (4.2.6), it follows by the chain rule that (_/)'(£)| is uniformly bounded in X n for each 
£. Therefore the arc length of each /)(£) is uniformly bounded above and below. Since there 
are 2 /i n such graphs, and 2 fi n ~ 2A n , the result follows. 

To prove (4.2.6), we note that by ( |4.2.5 ), 

2 7T 777. iy -1 ri 

_- 4 ( 0 ,£) 


\d e u n j(9,f)\ > C 2 X n e 


(0,0 £E X . 


We also note that from (3.2.4), d^Res^fl,£) = 0(0 where |£| < a(r) <C 1 in E Xn and 
^Ims^(04) = s'AO) + 0(£ 2 ) with s'A9 ) > Tj > 0. Thus, by another direct computation 


with the quasimodes (3.2.6), it follows that 

\d^u nJ (9,0\ < C 2 X n e 


h^-Im 4(0,0 


(0,0 £E K 


The analogue of (4.2.6) is therefore immediately true for the quasimodes u n j in place of 
the eigenfunctions u Xrt . Transferring to u Xri via Lemma |3.2.9| introduces errors that are 

2 7T 777. ^ j ^ 

0(e~ SXn / 6 ). However, in view of (4.2.5), these are lower-order than e L i 3 0 e _A ™ <5 / 10 
in E Xn . Such errors are therefore negligible, so the pointwise estimate in (4.2.6) is satisfied 
for the actual eigenfunctions u Xn . □ 

This gives both upper and lower bounds for nodal length near dominant boundary com¬ 
ponents, and the overall lower bound is an immediate corollary. It remains to estimate, from 
above, nodal lengths in the interior and near residual boundary components. 


4.3. Estimates for the harmonic extensions of interior eigenfunctions near resid¬ 
ual boundary components. To estimate nodal lengths all the way up to residual boundary 
components, we need to extend our domain slightly near the residual boundary components 
(independent of A) and use this slightly enlarged domain. In the following, we continue to 
work in the conformal model. So let D be an extension of D to include a (A-independent) 
neighborhood of the residual boundaries. We claim the following: 

Lemma 4.3.1. Given 5 > 0 as above, there exists an open domain D(5) D D with smooth 
boundary containing open neighbourhoods of the residual boundary circles such that the 
Steklov eigenfunctions tp x all have a harmonic extension to D(S ) denoted by u x . Moreover, 
in the extended neighbourhoods near the residual components one has the following estimate: 

I \ u *\ Ic 1 (5(5 )yd) — Ee 6 

for some geometric constant C > 0 depending on 5 but independent of X. 
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Proof. This follows from Proposition 2.4.1 Given any residual boundary component dDj 
and ho < consider the annular neighbourhood A,-(ho) D dDj with 

IIuaIIcHa^o)) < C , A 2 e A5 °||< / 2A||L2(aD i ) = 0(A 2 e A(_ 3 +(S o))). 

Choose h 0 < | and let D = D U (U^ =P+1 A,-(h 0 )). □ 

To simplify notation, in the following we denote the extended domain D(5) simply by D. 


4.4. Nodal length bounds away from dominant boundary components. We begin 
the analysis away from the dominant boundary components by considering simply connected 
sets in the interior of the extended domain D. In this section, we prove nodal bounds on 
these sets, then complete the proof with a covering argument. 


Lemma 4.4.1. Let U d U be simply connected open sets in D, each with C°° boundary, 
where: 

(1) U is not contained in Uf =l Aj(a), where a = a(r) was fixed in section f.2. 


(2) For each j — 1,. .., P, dU D dDj is a nontrivial circular subarc of dDj. 

(3) (Technical assumption) Translating the angle coordinate 9 if necessary, we assume 
that dU D Aj(a), for each j, may be parametrized in the form {(#,£); \9\ < 9 0 ,f = 
F(9)} with F(6 ) = 0 for \6\ < e 0 , F'(6) > 0 for 9 > e 0 and F'{9 ) < 0 for 9 < —e 0 . 

Then there exists a constant C, depending on U and U but independent of X, such that 

C(u\, U) < CX. (4.4.2) 


Remark 4.4.3. The technical assumption, made for convenience (see (4.4.8)), says that the 
boundary of U approaches each dominant boundary component dDj in a monotone fashion, 
travels along it for some distance, then departs in a monotone fashion. Note that F(-) 
measures distance from the boundary dDj. See Figure [l] for an illustration of sets U and U 
satisfying the hypotheses of the Lemma. 


Proof. Step 1: Pick xq G U \ Uj’ =1 Aj(a), which is possible by condition (1). Conformally 
mapping U to a disk and Xo to the origin keeps the functions harmonic, and changes the 
length of the nodal set by at most a geometric constant. Therefore, we may assume that U 
is a disk of radius 1 and that Xq is the origin. 

Note that U is contained in some disk B r with r < 1. We cover B r by finitely many disks 
Ai,... ,A k with centers Xi,... ,x k and radii ry,..., r k , where each Aj has the property that 
4A p the disk with center ay and radius 4ry, is contained in B\. By finiteness, we only need 
to prove the bound for each A,. By the nodal measure bound of Han and Lin 311. . Theorem 
2.3.1], 

2r- f |V«aI 2 

C(u x ,Ai) < CiN Ux (2Ai), where N^filAf) = 1 2Ai -^—. (4.4.4) 

Jd(2Ai) U \ 

Note that Han-Lin’s result is stated for a ball of radius 1, but scaling shows that it also 
holds for a ball of radius r, modulo a geometric constant which we absorb into C\. Since 
4Aj C Bi = U for all i = 1,..., k, a uniform control estimate on the frequency function 
([HX. Theorem 2.2.8], see also |NVl Section 3.2.2]) yields: 

N U x (2A*) < C 2 N Ux (U) 
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Then, summing up over the disks A; and using Green’s identity as well as the Cauchy-Schwarz 
inequality we get: 


C{u x ,U) < 


lu I^ u a| 2 | |«9zv^a | lat/1 11 dU 

sr^ £ 


8U U l 


IMG* ' 

lid := II ' llz, 2 (d)- 


(4.4.5) 


Here and further on we use the simplified notation 

Step 2: Since the argument for each dominant circle is the same, consider here the inter¬ 
section of the open set U with the r-annular neighbourhood of one fixed dominant boundary 
circle, say dD\. Let us estimate the numerator from above and the denominator from below, 
beginning with the denominator. We continue to divide through by a constant and without 
loss of generality assume that b\ l + + , _ = 1 , where &n,i,+ and b Hj i_ are constants defined 

Consider the norm over the portion of U consisting of an arc 7 along dD]_. 


3.1.1 


in Lemma 

By the interior quasimode approximation in Lemma 3.2.9, property (2) in Lemma 4.4.1 and 
a direct computation with the explicit quasimodes in (3.2.6), there is a constant c' > 0 such 
that 


IKHjU > 


I - 112 

1^71,1 | | ■y 


—c A 


(4.4.6) 


where u n ,i(x) is the global interior quasimode defined in (3.2.6). Let the curve segment 


7 = {(#,£)) l^l < £ o,£ = 0} with some e 0 > 0. Since Irn sf(9, £ = 0) = 0, it follows from 
Lemma 13.2.91 that 


|lln,l I | ry 


'| 6»|<£ 0 


27T 27T 

~rbn,i,+ cos(m n ,is( 0 )) + — 6 „, 1 - sin(m nj i s(9)) 

Li Li 


d0 > C(e 0 ) > 0 . 


Together with (4.4.6) this yields a lower bound 

II^aII^ > C'(£ 0 ) > 0 . 


(4.4.7) 


To get upper bounds for the numerator in (4.4.5), we split the integral over dU into three 
pieces: dU = 71 U 72 U 73 . The first piece, 71 = dU D A(a) is the putative leading term, 
coming from the part of dU inside the annulus A(a) around the dominant circle 8D\. As a 
consequence of our technical assumption (3), 


I Kill = 


'|0|<eo 


2tt 2tt 

~r bn,i,+ cos(m n ,is( 0 )) + — sin(m nj i s(9)) 

Li Li 


d6 


Go>| 0 |>eo 


O(e -r 1 A[F W+ O(F(0)3) ]) dQ + 0(e -AJ/2 } ^ 1? 


(4.4.8) 


since J go> \ 0 \ >£Q 0(e~ riX ^ e ' )+0< ' F ^ e ' , ' i ^) d6 = 0(A _1 ) for S 0 > 0 small by a change of variables 
9 e-)- F(9), using that \F'(6)\ > 0 in the range e 0 < |0| < k 

The second piece 72 = dU fl [D \ A(ct)]. This is the piece of dU in the closed domain D 
outside the annulus. Here, we know by formula (4.2.2) that with a = Ci(t) > 0 sufficiently 
small, 

||wa || 7 2 = ll«n,i|| 72 + 0(e 6 ) = 0(e G 


0(e~ 


Here, the first term is bounded using (3.2.6), (3.2.4) and (3.2.5) 
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Finally, the third piece 73 = <927 \ D is the p art of the boundary of 27 that is the exterior 
to D in the extended domain D. Using Lemma 


4.3.1 


we have 


IMF = 0(e- sx, ‘). 

Consequently, 

II^Allap = 0(1). 

For the normal derivative term, the same decomposition argument shows that 

I \9uU\\\ a ij = 0(A). 

Then, from (4.4.7), ( 4.4.9| ) and (4.4.10) it follows that 


£feir)<C l|4 “! lk !“ ,k < OA. 


IKII 2 


dU 


This completes the proof of the lemma. 


(4.4.9) 

(4.4.10) 

(4.4.11) 
□ 


Now we use a covering lemma to finish the proof. 


Lemma 4.4.12. Suppose LI is a connected, compact, smooth Riemannian manifold with 
smooth boundary <9f2 consisting of k connected components M x ,..., M k . Then there exists 
a covering of 0 by finitely many open s ets, ea ch of which is connected, contractible, and 


satisfies the assumptions on 27 in Lemma 4-4-1 


Proof. First take a finite collection of open balls (and half-balls, near the boundary) {27*} 
which covers k2 (compactness guarantees finiteness is possible). Take small enough balls to 
ensure that 27* is contractible for each i (for example, assume that each ball has radius less 
than the injectivity radius of 12). For each i, we will find an open set V) which contains 27*, 
has contractible closure, and moreover intersects each boundary component. The collection 
{Vi} is then the desired cover. 

So: begin with a set 27*. Since 12 is connected, for each j between 1 and k, let 7 j : [0, 1] — > 12 
be a smooth path with 7 ^( 0 ) G 27* and 7 ^( 1 ) G Mj. 

First, let t\ = supjt G [ 0 , 1 ]I71 (t) G 27 *}. Then let 71 be the portion of the image of 71 
corresponding to t G [t\, 1 ]; let W\ be open such that W\ is a contractible closed neighborhood 
of 71 , and let 27*. 1 = Ui U W\. Now 27*.i is contractible (as W\ contracts to 71 , and U * U 71 
contracts to Ui), and intersects Mi. 

Then proceed inductively to define 27*j; for each j, let tj = supjt G [0, 1] | 7 y (7) G 27*j_i}, 
let Tj be the portion of the image of 7 j corresponding to t G [tj, 1], and let Wj be the 
interior of a contractible closed neighborhood of 7 j. Then let 27* j = f / 7 ;,7 _ 1 U Wj. At each 
step, Uij is contractible, as it contracts to Uij -1 which is assumed contractible by the 
inductive hypothesis. And Uij intersects Mj. Appropriately smoothing out and adjusting the 
boundary of 27^ in each neighborhood Aj(a) gives us a new set V) satisfying all assumptions 
of Lemma |4.4.1 completing the proof. Note in particular that the technical assumption (3) 
is essentially local and can always be arranged. □ 


in Lemma 


4.4.1 


Applying this lemma to D, we obtain a finite collection of open sets 27*, i = 1, 
which cover D. For each x G D \ U ? =1 Aj{a), which is an open subset compactly contained 
in the interior of D, there is an open ball U x about x compactly contained in some Ui. 
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By compactness, a finite collection of these balls, {Uk} l k=1 , cover D \ U ? =1 Aj(a). Applying 
to each such Uj~ and the U t it is contained in gives us, for each k, 

£{Z UA fl Uk) < CkX. 


Lemma 


4.4.1 


Finiteness of the covering completes the proof of Theorem 1.2.1 


4.5. Real analytic Riemannian surfaces: general case. We consider here the general 
case when (Q, g ) is a compact real-analytic Riemannian surface with boundary dfl, and show 
that the argument above applies to this setting. It is well known that given a real-analytic 
Riemannian manifold it is always possible to extend it across its boundary to a slightly larger 
open real-analytic Riemannian manifold (see, for instance, Ed Example 5.50]). Let (12, g ) be 
such an extension of the surface (0, g). On the open Riemann surface 0, by the Behnke-Stein 
Theorem ([V] Theorem 9.4.6), there exists a strictly-plurisubharmonic exhaustion function 
for 0, and so there is a Green’s function G(z,w) G 7T(f2 x 0) with the property that 

A g>z G(z,w) = 8(w — z), ( z,w ) efixSl, 


and 


G(z,w) = --^log 


z — w\ — a(z, w), 


(4.5.1) 


where a is harmonic in both variables and hence a G C U (Q x 0). 

Next, we note that although we have used the Koebe conformal model in the genus zero 
planar case for convenience, it is not necessary. Indeed, returning to the initial Steklov eigen¬ 
functions ip\ G C u {dFt) one can directly derive the potential layer equations used to obtain 
the crucial quasimode approximations for the Steklov eigenfunctions along dFt. Therefore, 
the arguments of section A3 may be repeated with the only difference that one needs to take 
into account the contribution of a(z,w). However, since this function is harmonic, its contri¬ 
bution can be absorbed in the last two terms of the expression in Claim 2.3.7 Consequently, 
just as in the planar case, 

2 ^A n = X n Gcp Xn - \ n Kiip\ n - K 2 (p Xn , 

where the Schwartz kernels of K\ and K 2 are elements of C UJ (dfl x diT) and G is the same 


as in Claim 2.3.7 The quasimode approximation results in Proposition 2.3.4 and Lemma 


3.1.1 and the estimates in Lemmas 3.2.9 and 4.3.1 then follow as in the multiply-connected 


planar case after decomposing the boundary dQ into dominant and residual components. 
This completes the proof in the general case. □ 


Appendix A. Proofs of Lemma 13.1.11 and Proposition 13.1.31 


Throughout the appendix, for simplicity, we work with a real basis of eigenfunctions. 


A.l. An auxiliary lemma. Lemma 3.1.1 depends on the presence of recurring spectral 
gaps, so in order to prove it, we must first prove the following lemma. 


Lemma A. 1.1. Fix any sufficiently small e > 0. Then there exist pairwise disjoint closed 
intervals = [Aj, Bf\ C M + , i — 1,..., oo, with the following properties: 

(1) cr(P) U a(Mj) C \J*L x Ti; 

(2) A i+ i > Bi + e for all i; 
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(3) Fori > 2 and all j, each It contains at most one distinct element of cr(Mj) 

(4) For each i, there exists n such that X n G X \ and /i n G X\. 


Proof of Lemma \ A. 1.1 . Pick £ < j;, where L is the maximum of the boundary lengths 

Li, ..., L k . Since |A„ — fi n \ —> 0, there exists N such that if n > N, then \ n — fi n < e/2. 

Observe that since {/i n } is a union of k arithmetic progressions, each with period > 2i t/L, 
there are at most k elements of {/i n } in any interval of length 7 t/L. Thus any such interval 
must contain a gap of length at least 7 \/kL > 2e with no elements of {tin}- We may therefore 
choose some m 2 so that m2 > N and /i m2 > Pm 2 - i+2e. Consider the interval [/x m2 , /i m2 +7r/X] 
and observe that it must itself contain a gap of length at least 2e. Let n 2 > m2 be such that 
/i n2 is the left endpoint of the first such gap; we have fi n2 — Hm 2 < n/L. Then let m3 = n 2 + l, 
so that /i m3 is the right endpoint of that gap. Choosing n 3 as with ri 2 and iterating this 
process, we produce m 2 , n 2 , m 3 , n 3 ,.... 

Now let 

X\ = [0, Hm 2 ~ 3e/2]; X, = [/r mi - e/2,^ ni + e/2] Vj > 2. 

We claim that these intervals satisfy each property we want. Indeed, by construction, prop¬ 
erty (2) is automatic. Property (1) follows immediately from the fact that |A n — fx n \ < e/2 
whenever n > m 2 — 1 > N. To see property (3), note that each X t with i > 2 has length at 
most 7T /L + e < 27 t/L, and hence contains at most one element of each arithmetic progression 
with period > 27 t/L. Finally, property (4) is also immediate by construction (note that X\ 
contains /Xo — 0). This completes the proof. □ 

We now use these intervals to split the sequence {A n } into pieces with gaps of size at least 
e between each. Specifically, fix some e and let Xj, i — 1, 2 ,... be the intervals constructed 
int Lemma [A.1.1 For each i, we let Ci be the set of all j for which A j G Xj, and we say that 

j\ ~ j 2 4=^ 3i G N s.t. j\ G Xj and j 2 G Xj. (A.1.2) 

Since each interval Xj for i > 2 contains at most k eigenvalues X n (and each interval, including 
Xj, contains at least one), there exists a universal constant C such that 

j G Xj => j < Ci and i < j. (A. 1.3) 

A.2. Completing the proofs. 


Proof of Lemma \3. 1.1 . Let {e n } be an orthonormal basis of eigenfunctions for M, with eigen¬ 
values /i n , and note that each e n j = e n \ qd is either a trigonometric polynomial with fre¬ 
quency n n or identically zero. Write aj tk = (px-p Xfc)- By orthonormality and completeness 
of the eigenbases, we have 


^ dj,kP\j 
3 = 1 


Xa k 


a k,j e j- 


7=1 


We claim that if we set 


V'a„ — y ] a n,j€j, f An — 


a n,j&ji 


(A.2.1) 


(A.2.2) 




joon 


where ~ is defined by (A. 1.2), then the conditions of the Lemma are satisfied. Indeed, 
condition (1) is obvious, since by the definition of our intervals Xj, the frequency of each h,- 
with j ~ n is within 27 t/L of A n . It remains to prove (2). 
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In what follows, we let C and r > 0 be universal constants (which may depend on the 
geometry of D and on the conformal map from D to fl) and re-label at will. From Proposition 
2.3.4 and the Weyl asymptotics, we know that there exists r > 0 such that | Pe 3 - — < 

Ce~ rj . Plugging in the first equation in (A.2.1), we see that 


^ ^ 11 -J /Ij)^Anl|oO A 


~TJ 


n=l 

and hence that the same is the true of the L 2 -norm. Summing only over the n with n ^ j, 
we obtain 

£4(A-ft)’) l/ ’<Ce-’. (A.2.3) 

nooj 

By property (2) of Lemma A.l.l| we know |A„ — fij\ > e whenever n ^ j, so 

' 2 <Ce~ Tj . 




nooj 


n,j — 


However, by (A.2.1), for each j we have 'Yh n a ‘n,j = 1- So 


Ce~ Tj < J2 


< !• 


(A.2.4) 


(A.2.5) 




In addition, suppose j\ ~ j 2 with j\ < j 2 . Then, again using (A.2.1) and the orthonormality 
of <pk, as well as Cauchy-Schwarz, 

' 1/2 'E 


^ 1 a n,jx a n,j 2 I 


E 


a n,j 1 a n,j 2 I — 


<E 


n,j i/ 


n,j2 


) 1/2 < Ce~ Tjl . 


(A.2.6) 


n*j i nfj 1 n^ji 

Thinking of A = {a nj } as an infinite matrix, let Mi be the square submatrix of A with n 
and j in C . L ; note that by condition (3) of Lemma A.1.1 M t is of size at most 2k x 2k. We 
interpret (A.2.5) and ( A.2.6[ ) as saying that the columns of Mi have almost unit length and 
are almost orthogonal to each other. In particular, for each i, we may write Mf Mi = I + Ri. 
By (A.2.5) and (A.2.6), there exists a constant r > 0 such that 11-R*11oo < Ce~ Tl (the norm 


IlMHoo denotes the supremum of the entries of M). Therefore, for large enough i , (/ + R, L ) 
is always invertible, with inverse of the form / + Si, with HSiHoo < Ce~ Tt as well. In fact, 
(/ + Ri)~ l Mf Mi = /, so Mj -1 exists and equals (/ + Ri)^ l Mf = Mf + SiMf . Multiply the 
equation Mj Mi = I + R t on the left by Mj and on the right by M” 1 to yield 

MjMf = Mi(I + R^M - 1 = I + MiRi(I + SJM?. 

We see now that MiMj has the form / + Qi , where Qi satisfies ||<3i||oo < Ce~ Tl (note that 
the entries of Mj and Mf are bounded by 1). 

This now tells us that the rows of Mj have almost unit length: i.e. that 


1 - Ce~ Ti < 

j~n,jeCi 


<0 < !■ 


By subtraction and (A.1.3), there exists r > 0 such that 


E a h ^ Ce 

jo^n 


(A.2.7) 
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which is the analogue of (A.2.4), but summing in j instead of in n. Additionally, it is an 
immediate consequence of (A.2.4) and (A.2.7) that if j n, 

a n,j < Ce~ Tj ] al, < Ce ~™. 


n,j 


(A.2.8) 


By the Sobolev embedding theorem, for each k G N 0 , we have || • || c k {8D) < C\ \ • | \H k + 1 {dD)- 
We apply this to f\ n , knowing that we can compute Sobolev norms of the e 3 directly: 


II/. 


An | | (Jk 


< C\\f, 


n | IffH-i 2s 


< C ^(1 + l-lj) 


k +11 


n,j | 


J'yon 


By Weyl asymptotics of the /i j, 

ll/A.llc‘<C^(l+i )' !+1 " 2 - 


a n,j ~ 


C Y, (!+i)‘ +1 <i+C J] (1+j) 


k+1 a 2 ■ 

Ui n i j * 


joon 


j°°n,j<n 


joon, j>n 


Using (A.2.8), then the Weyl asymptotics again, we have 

+ C (1 +j) k+1 e~ Tj 


H/aJIc* < c ^2 (i + 


n )k+l e -rn 


j'#n,j<n 




< Cn( 1 + n) k+1 e~ Tn + C I (1 + t) k+i e~ Tt dt < CX L n +2 e 


k +1 —rt 


l k+2 —TTL 


(A.2.9) 


Choosing r slightly less than the current r/2, we can absorb and take square roots, 
completing the proof. □ 


Proof of Proposition 3.1.3 . From the discussion before (A.2.3), considering each term indi¬ 
vidually, we see that for every n and j, 

\a n j(X n . — Hj)\ <Ce~ Tj . (A.2.10) 

Using the matrix notation from the previous proof, let V t and W % be the diagonal matrices 
whose entries are {Aj : j G 22,} and {jij : j G respectively. Using (A.1.3), (A.2.10) 
implies the statement 

\\M i V i -W l M i \\ 00 <Ce- Ti . 


For sufficiently large i, Mi is invertible with inverse uniformly bounded in the 
norm. We deduce that 


matrix 


MiViMi = Wi + Ri, ||i2i||oo < Ce- 


So Wi + Ri is diagonalized by M t and has eigenvalues {Aj : j G £,}. By the Bauer-Fike 
theorem (see, e.g., [HJ, Observation 6.3.1]), the eigenvalues of Wi = (W* + Rf) — Ri lie in 


disks centered at each Aj of radius 


-i. 

OO 


\Ri 


loo < Ce~ Tl . 


Although the perturbed 
eigenvalues may move from one disk to another if there is overlap, there are at most k disks, 
so they can move at most 2 kCe~ Tl . Relabeling C = 2 kC, this shows that X n — /i n < Ce~ Tl , 
where X n G The result follows immediately by another use of (A.1.3), replacing r with 
n = t/C. ^ □ 
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